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PREFACE 


Pot Om thie contribution Of the Ontarl1o Ministry, or 
the Environment to the International Field Year for the 
Great Lakes (IFYGL) program was the estimation of ground- 
water inflow to Lake Ontario from the Canadian side, by 
extrapolating data from selected areas representative of 
larger hydrogeologic regions, This report, which describes 
the hydrogeology of the IFYGL Moira River, Wilton Creek and 
Thousand Islands study areas iS one in a series of reports 
dealing with the ground-water regimes of seven such 
selected representative areas along the Canadian shore of 
Lake Ontario, 


D,N, Jeffs, Director 
Water Resources Branch 


Toronto, January , 1981 
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ABSTRACT 


The assessment of the ground-water contribution to 
Lake Ontario from the hydrogeologic region along the north 
shore of the Lake, extending eastward from the Trent River to 
the eastern boundary of the Lake Ontario drainage basin on the 
north shore of the St. Lawrence River, was undertaken by the 
Ontario Ministry of the Environment as part of its contribut- 
ion to the International Field Year for the Great Lakes 
(IFYGL) program. Data were extrapolated from the hydrogeologic 
evaluation of the IFYGL Moira River, Wilton Creek and 
thousand Islands study areas which were considered to be re- 
presentative of the ground-water regime in this hydrogeologic 
region, 


Field investigations of the geology, data from water- 
well records on file with the Ontario Ministry of the Environ- 
ment and published literature were utilized for the assessment 
of the hydrogeology in the study areas. 


Ground-water supplies are obtained from the upper part 
of the bedrock and from the overlying Wisconsinan drift. The 
bedrock consists of crystalline rocks of Precambrian age over- 
lain by clastic rocks of Cambrian and Middle Ordovician age 
which are in turn overlain by marine carbonates of Middle 
Ordovician age. Well-productivity data suggest that the yield 
of bedrock wells is independent of the lithology and is 
primarily dependent on the fracture density in the bedrock. 
Water wells completed in the overburden are more productive 
than wells completed in the bedrock. 


An overall decrease in the storage of the ground-water 
reservoir occurred in the hydrogeologic region during the 
FreldsvYear period, (April J5.1972- to March 31, 1973),.as.identi- 
fied by the decline of water levels ranging from 0.18 to 0.97 
feet. This decrease in storage occurred although an increase 
iMeprecip1tation,.ranging -irom *3.07¢to 6.47, inches, was 
recorded over the long-term, annual average for the study areas. 


liyGraulic conductivity values, derived from transmissivity 
estimates obtained from short-term pumping tests, were assianed 
to the various materials in the study areas, From these values 
it was estimated that a grand total of about 44 cubic feet per 
secona of around water is discharging directly into Lake Ontario 
along the shore of this hy@e@rogeolodic region, 
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INTRODUCTION 


PURPOSE AND SCOPE 


Under the sponsorship of the Canadian and U.S. National 
committees for the International Hydrological Decade (IHD) 
program, the International Field Year for the Great Lakes (IFYGL) 
program was established to study the various hydrologic aspects 
of Lake Ontario and its drainage basin. Part of the Ontario 
Ministry of the Environment's contribution to the IFYGL program 
was the study of ground-water inflow to Lake Ontario, by extra- 
polating data from selected areas representative of larger 
hydrogeologic regions. 


The ground-water regime developed in any region is a result 
of the geology, topography, drainage and climate of that area. 
The eastern part of the Lake Ontario drainage basin, from the 
Trent River to the eastern boundary of the basin on the 
St. Lawrence River, is considered to be a hydrogeologic region 
where ground-water conditions are generally uniform. In this 
region, proceeding eastward from the Trent River, the Moira River, 
Wilton Creek and Thousand Islands study areas are situated on 
the north shore of the lake. 


For the purposes of the IFYGL study, to determine the amount 
of ground-water inflow to Lake Ontario, only the area adjacent 
to the lake was studied. Field investigations were made of the 
geology and test drilling was undertaken to provide information 
on aquifer characteristics, subsurface geology and ground-water 
levels in the area. Other aspects of the hydrogeologic assess- 
ment utilized information obtained from published literature and 
water-well records on file with the Ministry. 


LOCATION 


The Moira River IFYGL study area is located along the north 
and south shores of the Bay of Quinte in the counties of 
Hastings and Prince Edward; and extends eastward from the Trent 
River for approximately 17 miles to the Salmon River (Figure l). 
The study area extends northward for approximately three and 
one-half miles inland from the Bay of Quinte in the County of 
Hastings and approximately one-half mile south of the Bay in 
the County of Prince Edward. The Moira River enters the study 
area from the north at approximately latitude 44> 13' N. and 
longitude 77- 23' W. and then flows southward, draining into 
the Bay of Quinte at latitude 44° 09' N. 


The Wilton Creek IFYGL study area is located on the north 
shore of Lake Ontario, 14 miles west of the City of Kingston, 
in the County of Lennox and Addington (Figure 1). This study 
area is bounded on the south by Lake Ontario between longitudes 
76° 45' w. and 76° 53! W., a distance of approximately six miles, 
and extends inland (northwest) from the Lake for a distance of 
approximately six miles. Wilton Creek enters the study area 
from the northeast at latitude 44° 15' N. and longitude 76~ 49' 
W. and flows southwestward across the study area to drain into 
Hay Bay at latitude 44° 11' N. and longitude 76° 56' W. 


The Thousand Islands IFYGL study area is located along the 
north shoreuof the St, Lawrence River, 10 miles easty of tie 
City of Brockville in the United Counties of Leeds and Grenville. 
The study area extends approximately two miles inland from the 
St. Lawrence River and for a distance of 12 miles along the shore 
between Bucks Bay and Patterson Bay at longitudes 75° 59' W. and 
75° 49! W., respectively. 
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GEOGRAPHY 


PHYS IOGRAPHY 


Parts of three physiographic regions, the Napanee Plain, 
the Prince Edward Peninsula and the Leeds Knobs and Flats, as 
delineated by Chapman & Putnam (1966), are found in the three 
study areas (Map 1). The Moira River study area is comprised 
of the Napanee Plain on the north side of the Bay of Quinte and 
the Prince Edward Peninsula on the south side of the bay. The 
Napanee Plain comprises all of the Wilton Creek study area and 
the Thousand Islands study area is wholly in the Leeds Knobs 
and Flats physiographic region. 


Napanee Plain 


According to Chapman & Putnam (1966), the Napanee Plain is 
a flat to undulating plain of limestone which was scoured, 
Scraped and stripped of any existing overburden by the glaciers 
during Pleistocene time. The soil is generally only a few inches 
deep over much of the physiographic region. Some thicker deposits 
efsglacital till occur in the Moira River basin, including a few 
scattered, long and thin drumlins. The low areas in both the 
Moira River and Wilton Creek study areas are often found to be 
covered with a shallow veneer of stratified clay (Chapman and 
Putnam, 1966). 


Prince Edward Peninsula 


The Prince Edward Peninsula is the extension of the 
Napanee Plain across the Bay of Quinte, projecting into the 
eastern part of Lake Ontario. Similarly, low areas within the 
study area are mantled by a thin cover of stratified clay. 


Leeds Knobs and Flats 


The Thousand Islands study area is wholly enclosed by this 
physiographic region and consists of knobs of Precambrian rocks 
separated by flat lands of glacial drift (Chapman & Putnam, 1966). 
The flat lands are the result of glacial scour along zones of 
least resistance in the bedrock where glacial debris and other 
sediments were subsequently deposited. The rock knobs are 
generally bare as a result of wave action from the Champlain Sea 
which removed some of the shallow drift cover during the close 
of glaciation on the North American continent. A thin veneer 
of clay, deposited by the Champlain Sea, is usually found 
levelling the low areas between the Precambrian outcrops. 


TOPOGRAPHIC GRADIENTS AND RELIEF 


In the Moira River study area the regional topographic 
gradient is in the order of 33 feet per mile (0.0065 feet per 
foot) on the north side of the Bay of Quinte. In the’ vicinie, 
of the south shore of the Bay of Quinte (County of Prince Edward) 
the average topographic gradient is in the order of 75 feet per 
mile (0.0142 feet per foot). Maximum relief is approximately 
150 feet as measured from lake level at 245 feet above sea level 
(asl) to high land at 400 feet asl near the northern and southem 
boundaries of the study area on either side of the Bay of Quinte. 


In the Wilton Creek study area the regional topographic 
gradient is in the order of 33 feet per mile (0.0065 feet per 
foot) towards Lake Ontario. In the northwestern part of the 
study area the average slope of the land is southwestward to 
Hay Bay at approximately 20 feet per mile (0.0037 feet per foot). 
Maximum relief is approximately 150 feet as measured from Lake 
Ontario at 245. feet “asl to the tops of rock’ ridges in the norte 
part of the study area at 400 feet asl. 


The regional topographic gradient in the Thousand Islands 
study area is in the order of 40 feet per mile (0.0075 feet per 
foot) towards the St. Lawrence River in the south east. Locally, 
topographic gradients in the vicinity of the St. Lawrence River 
are in the order of 200 feet per mile (0.0378 feet per foot) as 
a result of the bedrock-outcrop pattern in the study area. 


DRAINAGE 


Drainage of the study areas is a result of the physiography 
where the regional slope of the land is towards the Bay of Quinte, 
Lake Ontario and the St. Lawrence River. With the exception of 
Prince Edward Peninsula, the major streams in the hydrogeologic 
region drain southward, generally having their headwaters in the 
Précambrian Shield to the north. South of the Bay of Ouinte, 
most streams drain to the north from the limestone highlands in 
the central portion of the Prince Edward Peninsula. 


Moira River Study Area 


Although the Moira River study area is part of the physio- 
graphic region of the Napanee Plain, drainage patterms in this 
area are different from the rest of the Napanee Plain to the 
east. In the Napanee Plain the dominant southwest direction 
of streamflow results from preferred glacial scour of the 
fracture: or jointing pattern developed! in the bedrock. =Structura: 
control from the bedrock is no longer dominant in the 
Moira River study area as evidenced by the southward-flowing 
pattern of creeks including the Moira River. This change can 
probably be attributed to the accumulation of glacial debris 
which is relatively thin but sufficient to alter the drainage 
pattern in this part of the study area. The prominent south- 


west direction of drainage, that is typical of the rest of the 
Napanee Plain, is present in the physiographic region of the 
Prince Edward Peninsula on the south side of the Bay of Quinte, 


Wilton Creek Study Area 


Drainage in the Wilton Creek study area can be divided 
into two parts: streams draining southeastward and directly into 
Lake Ontario in the south and streams draining southwestward 
into Hay Bay in the north part of the study area, Water from 
Hay Bay reaches Lake Ontario by way of Adolphus Reach which 
separates the physiographic regions of the Prince Edward 
Peninsula and the Napanee Plain, 


The development of long and narrow drainage basins parallel 
to each other in the study area, as well as over most of the 
limestone terrain of the Napanee Plain, appears to be related to 
a fracture or jointing pattern developed in the bedrock, The 
prominent southwest direction of streamflow was probably enhanced 
by glacial scour resulting from glacial movement in that 
direction, The southwest-northeast trending shore of Lake Ontario 
in the study area is also probably a development of structural and 
glacial control, Direct drainage to the Lake along the Lake 
Ontario shore is approximately at right angles to the shore, 
suggesting that this drainage pattern (southeast) is a consequence 
of structural control in the underlying bedrock, The largest 
stream in the study area is Wilton Creek which flows southwestward 
across the northern part of the study area, parallel to the Lake 
Ontario shore, This creek, as well as all the others in the 
north part of the study area, empties into Hay Bay. 


Thousand Islands Study Area 


The Thousand Islands study area shows an example of Pre- 
cambrian drainage where present-day streams flow in the low areas 
between ridges and knobs of bedrock, The primary trend of streams 
and their tributaries in the area is northeast-southwest, 
parallel to the St, Lawrence River and the inferred direction of 
glacial movement, Changes in the primary trend of stream 
direction generally occur at right angles and continue for short 
distances until the primary trend is re-established, The main 
streams eventually drain into the St. Lawrence River at places 
where the land surface is sufficiently subdued to permit a change 
in stream direction to allow drainage into the St, Lawrence, 

For example, the largest stream draining through the study area 
is La Rue Creek which empties into the St. Lawrence River at 

La Rue Mills, Most of the tributaries to the creek are either 
northeast or southwest flowing, parallel to the St, Lawrence 
River. In several places, the downstream portion of a 
tributary is parallel to the upstream portico of that tributary but flows 
in an opposite direction. This reversal of drainage suggests that drainage 
pattems in the area are influenced to a large degree by structural features 
associated with the Precanbrian bedrock. 


CLIMATE 


Data from meteorologic stations (Map 1) in the vicinity of 
the study areas are assumed to indicate the climate of the 
respective areas. The normal monthly and annual precipitation 
and temperature values for three representative stations are 
shown in tables 1 and 2 below. These long-term data (minimum 
23 years of record) indicate that the climate of the study areas 
(i.e. the entire hydrogeologic region) is relatively uniform. 
The climate of the hydrogeologic region appears to be slightly 
wetter towards the east as indicated by long-term, annual 
precipitation values of 33.86 inches at the Belleville meteoro- 
logic station (Moira River study area), 35.38 inches at the 
Kingston, Ontario Hydro station (14 miles east of the Wilton 
Creek study area) and 38.09 inches at the Brockville station 
(10 miles east of the Thousand Island study area). The long- 
term annual temperature of the region is a uniform 45°F., as 
indicated by the records shown in Table2. 
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GEOLOGY 
BEDROCK GEOLOGY 
General Stratigraphy 


The Moira River and Wilton Creek study areas are underlain 
by Paleozoic rocks of Middle Ordovician age which have been 
described by Liberty (1960 and 1971). The rocks in the 
Thousand Islands study area are of Precambrian age and have 
been described by earlier workers, notably Wynne-Edwards (1962, 
1963, 1967) and Hewitt (1963). Surface exposure of these rocks 
are found in creek valleys, in bedrock plains near Lake Ontario 
and in the bluffs along the Lake Ontario shore, Table 3, after 
Hewitt (1963) and Liberty (1960a, 1960b and 1971), indicates 
the stratigraphic sequence and includes a brief description of 
the lithology of the formations present in the study areas, 


Precambrian Lithology 


Rocks of Precambrian age outcrop in the Thousand Islands 
study area and are present at depth in the Moira River and 
Wilton Creek study areas. Since their deposition, the Pre- 
cambrian rocks have been subjected to folding, faulting, meta- 
morphism and other tectonic processes which have deformed the 
Original sediments into a complex system of hard crystalline 
rocks. 


The oldest rocks are the quartzites and gneisses of the 
Grenville Series whose distribution is shown in Map 2, These 
rocks were intruded by massive, dark-coloured gabbro, diorite 
and their altered equivalents (Buckingham Series) which outcrop 
in the north part of the Thousand Islands study area (Table 3). 
The Buckingham Series was in turn intruded by a series of 
granites and their respective oneisses (Table 3) which were 
later cut by dikes of diabase or porphorytic andesite, The 
complex structural and stratigraphic arrangement of the Pre- 
cambrian rocks is beyond the scope of this report and further 
reference is made to investigations conducted by Wynne-Fdwards 
G96 278 1963 7and 219.67), 


In summary the complex history of the Precambrian sediments 
has resulted in extensive changes in their original character 
so that they now appear as a relatively-homogeneous series of 
hard , crystalline, and well-jointed rocks, For the purposes 
of this report, the Precambrian aged rocks are considered to 
represent a single stratigraphic unit, 


Paleozoic Lithology and Depositional Environment 


Overlying the Precambrian surface are Paleozoic sediments 
which form the bedrock surface in the Moira River and Wilton 
Creek study areas, The lithologies described in the composite 
stratigraphic sequence (Table 3) for the study areas suggests 
that during Cambrian and Early Ordovician times, the Precambrian 
Shield was inundated by shallow, westward transgressing seas 
in which sediments of sand and clay were deposited, These 
basinal deposits are a lithostratigraphic unit comprised of the 
Potsdam and Shadow Lake formations of the Basal Group 
(Liberty,1969) and are in turn overlain by marine carbonate 
rocks; indicative of deepening water conditions during Middle 
Ordovician time, 
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The carbonate sequence of rocks (Simcoe Group) has been 
subdivided by Liberty (1960a, 1960b, and 1971) into the follow- 
ing formations in ascending order: Gull River, Bobcaygeon, 
Verulam and Lindsay. 


The Simcoe Group, as defined by Liberty (1969), is a rock 
unit consisting of limestone with no significant natural break 
in the lithology of the carbonate sequence within that group. 
The term Bobcaygeon Formation was introduced by Liberty (1969) 
for a lithogenic unit that lies between the lithographic lime- 
stone of the Gull River Formation and the equally distinctive 
interbedded limestone and calcareous claystone of the Verulam 
Formation. The Lindsay Formation overlies the Verulam Formation 
and consists of the remaining limestone sequence of the Simcoe 
Group. The areal distribution of these Paleozoic sediments in 
the study areas is shown in Map 2. 


The marine carbonates that comprise the Gull River, 
Bobcaygeon, Verulam and Lindsay formations indicate deepening 
conditions of the Middle Ordovician sea. Liberty (1969) suggests 
that the lithographic limestones of the Gull River Formation 
were formed from precipitated lime-mud. Shallowing conditions 
are indicated by the clastic nature (calcarenites) of the Bob- 
caygeon Formation. Fluctuating conditions of shallow and deeper 
water are indicated by the alternating limestone and claystone 
units of the Verulam Formation. The increase in calcareous 
material of the Lindsay Formation over the underlying Verulam 
Formation suggests slightly deeper and more stable conditions 
(Liberty, 1969) although fluctuating conditions are indicated by 
changes in sediment composition and the presence of rubbly 
materials. 


Paleozoic Structure 


The Paleozoic strata exhibit a gentle regional inclination 
that varies from 15 to 22 feet per mile to the south in the 
Moira River and Wilton Creek study areas (Liberty, 1960a; 1960b 
and 1971). The present attitudes of the strata are believed 
(Liberty, 1969) to be essentially those of their initial depo- 
sitional dip. Higher dips from eight to 15 degrees occur and 
probably reflect the drape of the Paleozoic sediments over the 
underlying high areas (monadnocks) on the Precambrian surface 
CI DerLy.. Loy). 


The Frontenac Axis is a ridge or platform of Precambrian 
rocks extending southeastward through the Thousand Islands study 
area and across the St. Lawrence River into the United States. 
The Paleozoic-Precambrian unconformity can be readily observed 
along the southwest flank of the Frontenac Axis (between the 
Thousand Islands and Wilton Creek study areas) where flat-lying 
sandstones of the Potsdam Formation uncomformably overly the 
Precambrian deposits. Locally, the absence of both the Potsdam 
and Shadow Lake formations in the Kingston area suggests that a 
local, Precambrian basement high existec during the Cambrian 
andinto the Ordovician (1.6. non deposition) at this locality 


(iiberty, 1971). 


Ait 


The Peterborough Arch is a name proposed by Liberty (1969) 
for a broad Precambrian platform extending southeastward from 
the Lake Simcoe area through the Moira River study area into 
the County of Prince Edward. This broad feature separates the 
Cambrian sediments on the east side of the arch in the Kingston 
area (Potsdam Formation) from those on the west side in the 
Trenton to Toronto area. Liberty (1969) also indicates. that ie 
Middle Ordovician carbonate sequence thickens southward, east of 
the arch. Several cisconformities have been delineated by 
Liberty (1963, 1971) in these Middle Ordovician strata that are 
attributed to the presence of the arch at the time of deposition. 
These disconformities occur within the Gull River Formation, 
between the Cull River and Bobcaygeon formations and at the base 
of the Verulam Formation. The presence of this arch is 
considered to be the reason for the alternating conditions of 
deeper and shallow water in the Lindsay Formation (liberty) 967 je. 


Faulting of the Paleozoic strata has been recognized in the 
Study areas. Most faults exhibit small throws of two to ten 
feet (Liberty,1971) although normal faults with 100) fect, of cine. 
have been recorded at Picton (Liberty, 1960b) and along the 
Salmon River (Liberty, 1963) in the Moira River study area 
(Map 1). This faulting is considered to be very recent in age 
for the most part. 


Bedrock Topography 


The elevation of the bedrock surface for the study areaceis 
ivbustrated in Map 2. The attitude of thie Ssurtace was coteined 
from records of water wells, on file with the Ministry of the 
Environment (Map 3), that penetrate the overburden to the bedrock. 
These data indicate that the bedrock Surface 1s similar to that 
Orethe present-day topographic or sland surtacesineall es three 
Study areas. A veneer of glacial debris is found capping the 
bedrock over parts of the areas but for the most part, bedrock 
outcrops are found very close to the land surface. 


PLEISTOCENE GEOLOGY 


The unconsolidated materials overlying the bedrock in the 
Study areas were deposited during the Pleistocene or Glacial 
Epoch and are considered to be Wisconsinan in age. Investiga- 
tions of these deposits have been made by earlier workers, 
notably Coleman (1936), Gillespie et al (1962), Gillespie et al 
(1963), Chapman and Putnam (1966), Henderson (1970), and 
Boni 197i) 


These Pleistocene deposits are described as being a thin 
veneer of glacial, glaciofluvial and glaciolacustrine materials 
overlying the bedrock. The surface expression and areal extent 
of these deposits are shown in Map 4. 


Only the latest events of the last glacial stage 
(Wisconsinan) are represented by the deposits, found in the: BYGp 
study areas. This suggests that glacial erosion was probably 
dominant in the eastern part of the Lake Ontario basin during 
most of the Wisconsinan Stage. Ice advanced from the east, down 
the St. Lawrence River into the Lake Ontario basin scouring and 
modifying the bedrock in its path. Eventually the ice filled 
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the Ontario basin and spilled out on to the surrounding high 
land, merging with ice from lobes advancing from the north to 
form a vast continental ice sheet over Ontario. Several major 
glacial advances and subsequent retreats occurred during the 
Wisconsinan Stage which lasted for approximately 100,000 years. 
The major events that occurred in the Great Lakes region have 
been documented by Dreimanis and Karrow (1972). 


FOr ther moctopange glacial till vor rubble is found-capping 
the bedrock, These deposits are generally less than 10 feet 
thick and are comprised of ground-up fragments of bedrock that 
vary from clay to gravel size. The till is generally described 
as being a stony sand till. Drumlins in the Wilton Creek study 
area indicate a southwesterly direction for ice movement whereas 
drumlins in the Moira River study area are oriented east-west, 
Suggesting a westerly shift in the direction of ice movement. 


Older deposits of sand and gravel of limited areal extent 
have been identified (Funk, 1977) below the till in the vicinity 
OfeWwiulton Creek and appear to be’ gqlaciofluvial inv origin. In 
the Thousand Islands study area, ice-contact deposits of sand 
and gravel have been delineated by Henderson (1970) and their 
extent also appears to be limited for the most part. 


At the close of the Wisconsinan glaciation, the ice front 
retreated to the east and north uncovering southerm Ontario. 
The Ontario Lobe withdrew from the Lake Ontario basin and stood 
in the St. Lawrence Valley. Meltwaters were ponded around the 
periphery of the ice forming high-level glacial lakes. As the 
ice front retreated, lower outlets for the glacial lakes were 
uncovered and the levels of the lakes dropped accordingly. Wave 
action in the lakes modified the underlying till plain and 
glaciolacustrine deposits of clay filled the depressions in the 
till plain. In the Moira River study area, thin deposits of sand 
(Map 4) are found near the present-day shore of Lake Ontario in 
the easter part of the study area and probably represent a water 
level slightly higher than the present lake (Coleman, 1936). 


As the glacier receded, an arm of the sea, called the 
Champlain Sea, entered the St. Lawrence-Ottawa Lowland. The 
western limit of the sea is not precisely defined but is believed 
to have extended into the eastern part of the Lake Ontario basin. 
Chapman & Putnam (1966) suggest that the western shore of the 
Champlain Sea was on the rocky edge of the Precambrian Shield 
(Frontenac Axis) where conditions were not favourable for the 
formation of beaches. Deposits of marine clay and sand have 
been delineated by Henderson (1970) in the Thousand Islands 
study area. Some of the clay sediments in the Wilton Creek study 
area are believed to have been deposited by the Champlain Sea 
(Gillespie et al, 1963). 


ae 


HY DRO GEOLOGY 


GROUN D-WATER OCCURRENCE 


Ground water is an aspect of the hydrologic cycle which 
has resulted from the infiltration of water at the earth's 
surface. The water has moved through the unsaturated zone 
into the saturated zone and, under the influence of gravity, 
to a point Of Natural Or artiticial, discharge pinoduct boar, 
lake or well. The water stored in the saturated zone is called 
ground water and natural recharge to this system generally 
occurs either directly from precipitation infiltrating into 
the ground or indirectly from snowmelt. 


An aquifer is defined as a water-bearing horizon that will 
transmit significant amounts of water to a well. In the 
exploration for ground-water resources, the connotation of the 
wording “Significant amounts" further defines the term aquifer 
with respect to the usage of water for domestic, industrial or 
municipal supplies. An aquifer in the context of this report 
is one that will supply water for limited domestic use. In the 
study areas, aquifers are found both within the overburden 
deposits and in the bedrock. 


Ground water is stored in the interstices or openings or 
pore spaces comprising the material through which the ground 
Water moves. The capacity Of any material to hold or yield 
water is dependent on the porosity and permeability of the 
material (i.e. the total number of interstices or voids and 
the interconnection of these interstices, respectively). 
Thus, a high porosity material such as shale or clay will not 
yield large quantities of water due to absorption phenomena 
resulting from the molecular size and poor interconnection of 
the pore spaces. 


Within the overburden, the most productive materials are 
the granular deposits of sand and gravel. These permeable 
materials, when encountered in the subsurface, are utilized 
as sources of ground water. These deposits, however, do not 
appear to be abundantly distributed in the study areas, as 
indicated from water-well records on file with the MinilStry 7o. 
the Environment and the surficial geology shown in Map 4. 


Numerous bored and dug wells penetrating low water-yielding 
materials such as glacial deposits of till or glaciolacustrine 
silts and clays are also exploited to supply water for domestic 
needs. The yields from wells completed in the above-mentioned, 
low-yield materials, are essentially controlled by the storage 
Capacity of the large-diameter wells. Recent studies in the 
prairies (Meyboom et al,1966) suggest that the movement and 
storage of ground water in glacial tills and some clays is within 
Joints or cracks developed in these materials. 


. The bedrock comprised of the marine carbonates of the 
Simcoe Group (Ordovician) and the crystalline rocks of the 
Precambrian, is generally a poor aquifer in the study areas. 
Microscopic examination indicates that the bedrock porosity is 
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low resulting from the past history of these rocks (compaction, 
consolidation, metamorphism, etc.) since their deposition. The 
rocks of the Cambro-Ordovician Basal Group (Potsdam and Shadow 
Lake formations) are better aquifers as a result of their 
granular nature. Although the primary porosity of these granular 
materials has been restricted by subsequent cementation of the 
grains, yields of water wells completed in these formations 
(Sobanski,1970) are substantially better than those wells 
completed in the Ordovician limestones or the crystalline rocks 
of the Precambrian, 


Secondary porosities consisting of joints, fractures, 
cleavage planes and faults were subsequently developed in the 
bedrock, probably as a result of tectonic processes. Other 
processes, such as those associated with weathering have enlarged 
these planes (including the original bedding planes), creating 
secondary voids through which the ground water presently flows. 
These secondary openings tend to become tighter and fewer in 
number with increasing depth, 


GROUND-WATER LEVELS 
Change in Storage 


The change in storage of the ground-water reservoir can be 
identified by the measurement of rising and declining water levels 
in water wells or boreholes. During periods of high precipitation, 
storage in the ground-water reservoir is increased after the 
soil moisture deficiency (SMD) in the unsaturated zone is re- 
plenished. When the SMD is satisfied, the infiltrating water 
will percolate downwards into the saturated zone under the 
influence of gravity. Conversely, storage is decreased under 
conditions where the infiltrating water is intercepted before 
reaching the saturated zone while discharge from the saturated 
zone continues. These conditions occur during periods of low 
precipitation and/or excessive use of subsurface water by man or 
vegetation, 


An observation-well network of private and abandoned wells 
was eStablished for the IFYGL program in the Moira River, 
Wilton Creek and Thousand Islands study areas. The change in 
water levels was computed for the wells in these networks over 
the wretaqvear period, from April 1, 1972) to March 317 4197/3 
and the results are presented in Table 4. 


The data in Table 4 suggest that an overall decrease in 
the storage of the ground-water reservoir occurred in all the 
IFYGL study areas during the Field Year. This is indicated 
Dytameave rage decline of (0597, 0.18 sand 0-26 steelrof water 
levels measured in the observation-well networks for the 
Moira River, Wilton Creek and Thousand Islands study areas, 
respectively. 
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Climatological data, presented in Table 5, indicate that 
excess precipitation over the long-term average occurred in 
all the study areas during the Field Year, Long-term data 
are not available for the Morven IHD station and it was 
assumed that long-term data from the Kingston Ontario Hydro 
station (approximately 14 miles east of the Wilton Creek 
study area) were representative for the Wilton Creek study 
area, Other IFYGL studies to the west, in the Bowmanville- 
Newcastle area, Duffins, Oakville and Forty Mile creeks, 
indicate that excess precipitation over the long-term average 
contributes to a rise in ground-water levels, The decrease 
in ground-water levels that occurred during a period of excess 
precipitation in the Moira River, Wilton Creek and Thousand 
Islands study areas suggests that the present ground-water 
resources are being slowly depleted by the consumption of 
water from the existing wells in the study areas, It would 
appear that a vrolonagec decrease in precipitation from the 
long-term annual average will prohably result in a shortage 
of ground-water supplies in these areas, 


Observation-Well Hydrographs 


Observation-well hydrographs at site WC-8 (Figure 2), in 
the Wilton Creek study area for the years 1968 to 1973, 
indicate that the maximum fluctuation of ground-water levels 
over this period of time is approximately six feet, Long-term 
hydrographs are not available for the Moira River and 
Thousand Islands study areas but a comparison of observation- 
well hydrographs from these areas over the Field Year period 
Eeomenprrlewds, 1972 fo March Sip 1993 (Figure 3) indicates wa 
Similar response over that period, Maximum peaks and lows 
in the water-level stages occur at the same time in the study 
areas, The lowest depth of water level was recorded in the 
month of September with a peak in late December - early 
January and aminor low at the end of February (Figure 3). 
This annual cyclic phenomenon results from the change in the 
rate of recharge to the ground-water reservoir, 
Climatological data for the study areas indicate that the 
long-term monthly precipitation (Table 1) is relatively 
constant, ranging from 2% to 34 inches per month, During the 
growing season (May to October inclusive) a decline in water 
levels occurs, reflecting a soil-moisture deficiency (SMD) 
created in the unsaturated zone from the consumption of 
available water by vegetation and evaporation from the land 
surface, A corresponding rise in water levels occurs during 
the non-growing season (November to April inclusive) after 
the SMD in the unsaturated zone is replenished and recharge 
to the ground-water reservoir resumes, The low recorded in 
Pebruary 1S probably due to freezing of the land Surface, 
consequently inhibiting recharge to the ground-water reservoir 
while discharge from the ground-water reservoir continues. 
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Piezometer at 1 76 feet 


Piezometer at 29 feet 
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Figure 2. Observation-well hydrograph for site WC -8in the IFYGL Wilton Creek study area 
(after Funk, 1977). 
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As indicated earlier, the primary source of ground water 
in the study areas is the bedrock which is generally found 
OULCLOpDING, Ob very close to the land surface. “If Tt is 
assumed that recharge to the ground-water reservoir in the bedrock 
is constant, the magnitude of the ground-water level 
fluctuations can be used as a qualitative measure of the hy- 
draulic properties. Large fluctuations in ground-water levels 
would represent limited storage capacity in the bedrock and 
conversely small fluctuations would represent large storage 
Capacity. During the Field Year period an excess of pre- 
Chota CLOMyaGangiuyg srom oe 7O"tor6.47 snchesr over the long- 
termyanniual average (Table 5), was recorded as occurring in 
the three study areas. The response of the observation wells 
to this input was similar (Figure 3); however,the magnitude of the 
water-level fluctuations for-the Field Year period ranged,from 
ee joetonVOrSG feet’ This variability in) magnitude suggests: taal 
Seqontetrcant differences in the hydraulic properties (i:e. storace 
Capacity) ef the materials’ are present’ (Figure 3) (aeigk ites 
assumed that a constant rate of recharge occurred over tne stucy 
areas. 


GROUND-~--WATER MOVEMENT 
General 


The ground-water flow field is a three dimensional field 
thatroccurs” From the tovovof the saturated zone down’ to "ai depth 
where ground-water flow becomes impossible (i.e. where the 
material in the saturated zone becomes completely impermeable). 
The ground water at any point in this field can be expressed in 
terms of total head defined’by Hubbert (1953) as the sum of 
pressure head, elevation head and velocity head at that point; 
however, the velocity head is so small for ground-water con- 
siderations, it can be neglected. The static water level (the 
level to which water rises) in a well is a measure of the 
ground-water head at the point where the water was found. The 
total flow vector (Meyboom et al, 1966) of the three-dimensional, 
ground-water field can be resolved into two, two-dimensional 
herizontalvand vertical vectors. The "flow resultant” as the 
flow component of the ground water in the vertical plane and 
the "horizontal component" is the flow component of the ground- 
water in the horizontal plane (Meyboom et al, 1966). In 
either case, the component of flow that is depicted (flow 
resultant or horizontal component) is assumed to essentially 
represent the total flow vector of the ground-water flow field. 


Water-Level Data 


As indicated above, the static water levels in wells are 
a measure of the potential causing ground-water flow. Water- 
level data from water-well records on file with the Ministry 
Ofetne Environment suggest that the static water level of any 
well, completed in any water-bearing horizon, is usually found 
within a tew feet of the land surface, regardless of the» season 
oer year in which the well was drilled. An average value for 
depth to static water level was computed for 866, 205 and 282 
bedrock wells in the Moira River, Wilton Creek and Thousand 
Islands study areas, respectively. The average, static 
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water-level value for wells in each of the above-mentioned 
study areas was found to be 1] feet (standard deviation of93 
feet), 14 feet (standard deviation of 9 feet) and 14 feet 
(standard deviation of 10 feet), respectively. Although 

these wells were drilled in different seasons and completed 

in different water-bearing horizons over a span of approximately 
20 years, the average, static water level is within a few 

feet of the land surface. This feature can be observed in 
figures 2 and 3 where the maximum depth to the static water 
level occurred in September, 1972 during the Field Year period. 
The water-level depths were 13 feet in the Moira River study 
area, 21 feet in the Wilton Creek study area and 11 feet in 

the Thousand Islands study area. Anomalously deep water- 
levels can usually be traced to excessive withdrawals or 
consumption of ground-water or to deep wells completed in the 
less active or stagnant part of the ground-water reservoir, 

All of these data suggest that a hydraulic connection exists 
within the bedrock and that the general direction of ground- 
water flow can be determined from these water-level data. 


Potentiometric Surface 


The water-level data from water-well records in the study 
areas were used to compile a potentiometric surface map for 
ground water in the bedrock. A potentiometric map was not 
prepared for wells completed in the overburden as insufficient 
data were available. The hydraulic-head values (water-level 
elevations) for the overburden wells are shown separately on 
the map (Map 5) prepared for the bedrock. ‘The similarity in 
the two sets of data, for the overburden and bedrock, suggests 
that they are hydraulically connected. For purposes of this 
report, the overburden and the uppermost part of the bedrock 
are considered to constitute the active part of the ground- 
water reservoir. 


Lines of equal potential were interpolated between values 
of hydraulic head obtained from water-well records on file 
with the Ministry of the Environment (Map 5). By definition 
ground-water flow is at a direction normal to the lines of 
equal potential. Examination of Map 5 indicates that the 
potentiometric surface conforms to the topographic expression 
of the ground surface and may be considered to be a subdued 
replica of that surface. The ground-water divides appear to 
be generally coincident with the surface-water divides and by 
analogy, areas of recharge and discharge are topographic highs 
and lows, respectively. Thus, ground-water flow in a lateral 
direction is controlled by the topography; movement being 
from topographically high to topographically low areas. 


Hydraulic Gradient 


Assuming that the "horizontal component" of ground-water 
flow (as illustrated in Map 5) represents the total flow 
vector of the ground-water flow field, the slope of the 
potentiometric surface is the hydraulic gradient (I) under 
which ground-water movement occurs. Ground-water movement 
is dependent on the hydraulic gradient (I) according to 


Ze 


Darcy's Law, Q = KIA, where: 

Quantity of flow per unit time 

Hydraulic conductivity 

Hydraulic gradient 

Cross~sectional area through which movement 
takes place. 


PHR O 


The hydraulic gradient (I) is the head loss (in feet of water) per unit 
distance (foot) of flow path, resulting from the frictional resistance 
to flow within the voids of the material throuch vhich the ground water 
moves. 


As indicated earlier, the potentiometric surface (Map 5) is similar 
to the land surface and ground-water flow is from high to low potential 
in a direction normal to the lines of equal potential. The hydraulic 
gradient (I) is in the order of 0.0065 feet per foot in the Moira River 
and Wilton Creek study areas and 0.0075 feet per foot in the Thousand 
Islands study area (Map 5). 


AQUIFER CHARACTERISTICS 


The water-yielding properties of an aquifer are a function of pore 
size and their interconnection. For quantitative ground-water studies, 
an estimate of the aquifer characteristics is required. The most common 
method for determining the water-yielding properties of an aquifer is by 
means Of a pumping test. The yield of a well, as determined from a 
pumping test, may be used to provide information on the water-yielding 
properties of the materials encountered such as specific capacity. 

The specific capacity of a well is its yield, in Imperial gallons per 
minute per foot of drawdown (Igpm/ft), for a stated pumping period and 
rate. The specific capacity is numerically expressed as O/s, where: 


Q 
Ss 


pumping rate (Igpm) 
drawdown in feet 


oil 


In addition, pumping-test data may also be used to provide an 
estimate of the transmissivity and the hydraulic conductivity of the 
material in which the well is completed. The transmissivity (T) is 
defined as the rate of flow of water in imperial gallons per day (IGPD) 
threugh a vertical strip of the aquifer one foot in width and extending 
the full saturated thickness of the aquifer under a hydraulic gradient 
of one foot per foot at the prevailing temperature of water. 

The hydraulic conductivity (K) may be determined using the relation 
K = T where: 

m 
hydraulic conductivity (ft/day) 
transmissivity (IGPD/ft) 
contributing thickness of the aquifer (ft) 


y 
tou 


Specific Capacity 


Statistical analyses of specific--capacity data (derived from pumping 
tests) by means of specific-capacity frequency graph (Figure 4) provide 
a convenient method of camparing the water-yielding properties of the 
material in the study areas. 
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Figure 4. Specific- capacity frequency graphs for bedrock wellsinthe IFYGL Moira River, Wilton Creek 
and Thousand Islands study areas. 
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Insufficient data were available to compile frequency 
graphs on the overburden materials, The specific-capacity 
data for wells completed in the bedrock were tabulated in 
order of magnitude and the frequencies were computed by the 
Kimball (1946) method. Values of specific capacity were 
plotted against percent of wells on logarithmic probability 
paper, The slope of the line passing through the data points 
indicates the variability of the water-yielding properties of 
the material; steeper slopes signify greater variability than 
flatter slopes. Comparison of the frequency curves for the 
three study areas (Figure 4) indicate a similarity in slope 
Ci,eG. vatblabliity). In addition, the productivity of the 
wells is also similar with 50% of the bedrock wells completed 
in the study areas having an average specific capacity of 
approximately 0.20 Igpm/ft of drawdown (Figure 4), 


These specific-capacity data suggest that the yield of 
bedrock wells in the study areas is independent of the 
lithology, The productivity and variability of the bedrock 
in the study areas appears to be the same (Figure 4), although 
wells in the Moira River and Wilton Creek study areas are 
completed in marine carbonates and claystones of Middle 
Ordovician age and wells in the Thousand Islands study area 
are completed in hard crystalline rocks of Precambrian age, 
This further suggests that the well yields are primarily 
dependent on the fracture density of the bedrock which is a 
result of tectonic processes to which the rocks have been 
subjected since their deposition, Specific-capacity data are 
not available in the study areas to permit a comparison with 
wells completed in the sandstones and arkoses of the Potsdam 
and Shadow Lake formations (Basal Group). Data on file with 
the Ministry of the Fnvironment (Sobanski) ,1970) indicate 
that these formations are more productive (specific-capacity 
of 3.5 Igpm/ft) than those comprising the bedrock surface 
(Simcoe Group and Precambrian) in the study areas. 


Transmissivity (Modified Non-Fgquilibrium Well Formula) 


Information from the short-term pumping test data from 
the water-well records on file with the Ministry of the 
Environment were utilized to estimate the transmissivity 
(and by analogy, the hydraulic conductivity of the different 
materials). Transmissivities were obtained using the modified 
non-equilibrium well formula (Cooper and Jacob, 1964): 


ris 264.0 
As 
where: T= transmussavicy CIGPD/£Le) 
Q = pumping rate in Imperial gallons 
per minute (IGPM), 
As = change in drawdown per logi9 


cycle (feet) 
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The derivation of the formula is based on the following 
assumptions (Johnson 1966): 


1) the water-bearing formation is uniform in character 
and permeability in both horizontal and vertical 
directions, 


2) the formation has uniform thickness, 
3) the formation has infinite areal extent, 
4) the formation receives no recharge from any source, 


5) the pumped well penetrates and receives water from 
the full thickness of the water-bearing formation, 


6) the water removed from storage is discharged 
instantaneously with lowering the head. 


The above-mentioned assumptions are not met in totality; 
however, it is assumed that the transmissivity values obtained 
from the well-record data are within an order of magnitude of 
the absolute values. 


The results of the computations of hydraulic conductivity 


for the various materials in the study areas are presented in 
Table 6. 
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Transmissivity (Specific Capacity) 


The transmissivity of ar aquifer can he estimated from the 
specific capacity of production wells. The theoretical specific 
capacity of a well may be written as: 


ae T 
114.6 Wu; L4/B) 
where: 
2 (2g 
Wee 242 ae S 
ce 
“w/B — W 
7 Pn 


W(u,r /2 )= well function for leaky artesian aquifers 
e Hantush, 1966), 


Q = discharge of pumped well (Igum), 
Ss = drawdown (feet), 


Yr = nominal radius of well (in feet), 
T = transmissivity (Igpd/ft) 

S = coefficient of storage, 

t = pumping period in minutes, 


p' = coefficient of vertical.permeabil ity of 
confining beds (Igpd/ft’), 


m' = saturated thickness of confining bed (feet) 


Reliable data are not available on the coefficient of storage 

(S) in the study areas. The magnitude of S depends on the 
elasticity of the aquifer material and the fluid. In confined 
aquifers j_ S does,not show large variations, generally ranging 
from 10~~ to 10-" (Kruseman et al, 1970). In the study areas, 
aimost all the: well. records: on file with the Ministry of the 
Environment indicate confined conditions (i.e. having a 

static water level above that where the water was found). 
According to Meyer (Bentall, 1963), any changes in the 
coefficient of storage (S) correspond to only small chances 

in the transmissivity (T) and specific capacity (0/s); 
tiererore:, inaccuracy in estimating S 1s mot a serious jlamicing 
Factor. 


With regard to the vertical permeability of the confining 


bed (P'), no data are available on this parameter in the study 
areas. 


28 


In addition to being dependent on the hydraulic properties 
of the aquifer and confining beds, the specific capacity (0/s) 
will vary as a result of the following factors: 


i) as (the nominal radius of the well) - O/s increases 
with increasing Cys 1.6. directly proportLondl 2o 
HOG. te 
WwW a 
2.) t (the pumping period) - Q/s decreases with increasing 


time, 1.e€. inversely proportional to log t, 


3% s_ (the well loss) '- Q/s decreases with an increase in 
pumping rate in wells with a high Sue 

The relationship between the theoretical specific capacity 
and the transmissivity of an artesian aquifer is shown in 
Figure 5(after Csallany and Walton, 1963). It is assumed that 
the conditions depicted in Figure 5 are applicable to the Moira 
River, Wilton Creek and Thousand Islands study areas. 
Hydraulic conductivities were computed from specific-capacity 
data, using Figure 5, for artesian wells (confined conditions) 
in the study areas and the results are presented in Table 7. 
Accoraing «to, Walton (1962) and Walton et: al (1962),, thestrans- 
missivity that is computed from specific-capacity data does not 
vary Significantly between a pumping period of one to eight 
hours and a limited range of the nominal radius (r._) of the well. 
In order to minimize errors, only wells with a pumping period 
ranging from one to ten hours and with a well diameter of four 
to eight inches were used in the computations of hydraulic 
Conductivity from specific-capacity-data (Table 7). 


The data in Table 7 indicate that the hydraulic conductiv- 
ities derived from the specific-capacity data are larger by 
approximately a factor of four to six than those hydraulic 
conductivities derived from the modified, non-equilibrium well 
formula. Sufficient data on well losses (s_) and hydraulic 
properties of the confining beds are not available to delineate 
the hydraulic conductivity from the specific-capacity data to, any 
greater degree of accuracy. 


GROUND-WATER CONTRIBUTION TO LAKE ONTARIO 
tnteroduction 


mie slarcy formula, O.= TIL, was used; to compute the anounc 
of direct ground-water inflow to Lake Ontario, where: 


quantity of water in imperial gallons per day (IGPD), 
transmissivity (IGPD/ft) 

hydraulic gradient (feet per foot), 

width of cross section (feet). 
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In order thet values be assigned for the parameters listed 
previously, several assumptions on the ground-water flow field 
in the study areas were made. The assumptions, with respect to 
contributing area, depth of active flow, saturated thickness, 
transmissivity and hydraulic gradient are discussed below. 


Contributing Area 


As indicated in an earlier section, the area near the Lake 
Ontario shore is considered to be in a zone of general ground- 
water discharge. In order to facilitate the computationsor 
ground-water inflow to the Lake, it was assumec that the area 
within a distance of one mile from the shore was representative 
of the conditions in the zone of ground-water discharge. 


Depth of Active Flow 


The mean value of bedrock penetration for wells completed 
in the bedrock was found to be 32, 41 and 69 feet, respectively, 
for the Moira River, Wilton Creek and Thousand Islands study 
areas. From the above data it was assumed that significant 
ground-water flow in the bedrock would not extend below 50 feet 
of bedrock penetration in the Moira River and Wilton Creek study 
areas, and below 100 feet in the Thousand Islands study area. 
The larger figure for bedrock penetration in the Thousand Islands 
study area probably reflects the harder nature of the 
crystalling rocks and suggests that fractures extend to a greater 
depth in these rocks than fractures developed in a carbonate 
terrain. These arbitrary values of 50 and 100 feet of bedrock 
penetration in the study areas are considered to be the limiting 
boundary between the active and stagnant parts of the ground- 
water reservoir. 


Saturated Thickness 


The reported, static water level in the well records was 
used to determine the saturated thickness of the materials 
within a distance of one mile from the Lake Ontario shore in 
the study areas: In many wells a relatively low, static water 
fevel in conjunction with a thin overburden thickness indicated 
Chat the contribution of ground water from the overburden wae 
negligible. 


Transmissivity 


Hydraulic conductivity values, derived from short-term 
pumping tests (Table 6) using the modified, non-equilibrium well 
formula (Cooper and Jacob, 1964) were used to compute the trans- 
missivities of the materials within the zones of interest. 

The various materlals reported in all the logs of wells within 

a distance of one mile from the Lake Ontario shore were assioned 
a hydraulic conductivity which was then multiplied by the 
Saturated thickness to determine the overall transmissivity in 
the active part of the ground-water reservoir. For those wells 
enaing in overburden an additional transmissivity for 50 feet 

of bedrock in the Moira River and Wilton Creek study areas,and 
for 100 feet of bedrock in the Thousand Islands study area was 
used. For the computations of the bedrock wells, it was assumed 
that all the wells penetrated the full thickness of the active 
part of the ground-water reservoir. 


D2 


The total transmissivities for each well were summed and the 
results are presented in Table 8. It was further assumed that 
the transmissivity values of about 2400, 1500 and 3400 IGPD/ft 
(Table 8) are representative for the shore area under 
consideration in this hydrogeologic region. 


Hydraulic Gradient 


The average hydraulic gradient (I) was assumed to be 
0.0065 feet per foot in the Moira River and Wilton Creek study 
areas and 0.0075 feet per foot in the Thousand Islands study 
area (Map 5) as indicated in the earlier section on ground-water 
movement. 


Computation of Ground-Water Inflow 

The values for T and I, as derived above, were substituted 
in the Darcy formula 9 = TIL, to obtain a value of ground-water 
contribution to Lake Ontario, per mile length of shoreline as 
follows: 


Moira River 


Ora TEL 


Za Joe OU OD oe 5260 


83,638 IGPD per mile length (of shoreline) 


0.155 cubic feet per second (cfs) ver mile length 
Pe ptonmCereel 


Ove LG 


WA 00D 1X O50 


49,935 IGPD per mile length (of shoreline) 


0.093 cubic feet per second (cfs) per mile length 
Thousand Islands 


Oya LE 


II 


Bio Coes 0 0075x5290 


134,046 IGPD per mile length (of shoreline) 


II 


0n249 “cubic feet per second: (cis) per mise: lenge) 
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The transmissivity values, as derived above, were used 
to compute the total volume of ground water discharging to 
Lake Ontario from the hydrogeologic region in which the Moira 
River, Wilton Creek and Thousand Islands are considered to be 
representative. This region extends eastward from the Trent 
River to the eastern boundary of the Lake Ontario drainage 
basin on the St. Lawrence River at the Town of Prescott 
(Figure 1). The Moira River study area is assumed to be 
representative of the westerm portion of the hydrogeologic 
region which extends eastward for 31 miles along the Bay of Quinte 
from the Trent River to the Napanee River. The Wilton Creek 
study area is considered to be representative of the central 
part of the hydrogeologic region that extends for approximately 
59 miles of shoreline from the Napanee River, along the east 
shore of Long Reach, the north shores of Adlolphus Reach, the 
North channel and Lake Ontario, to the Town of Gananoque on 
the north shore of the St. Lawrence River. The Thousand Islands 
study area is representative of the easter portion of the 
hydrogeologic region which consists of approximately 46 miles 
of shoreline along the north shore of the St. Lawrence River 
from the Town of Gananoque to the Town of Prescott. Using the 
above-mentioned values, it was estimated that 4.81, 5.49 and 
11.45 cfs of ground-water are discharging into Lake Ontario 
from the west, central and eastern portions of the hydrogeologic 
region, respectively. These values give a total of 21.75 cfs 
of ground-water discharge from the hydrogeologic region. This 
Valuewot 2 ojoncrs-1s Simivaruto the value Gl “Loto i cfs stor 
roughly the same area as computed by Haefeli (1972). 


Additional Ground-Water Inflow 


An additional amount of direct ground-water discharge to 
Lake Ontario can be computed for the shoreline in the County 
O£ Prince Edward which forms part of the Moira River study area. 
A total of approximately 144 miles of shoreline bordering the 
Bay of Quinte and Lake Ontario are present in the county. If 
it is assumed that the value of 0.155 cfs per mile length of 
shoreline in the Moira River study area is applicable to the 
area, an additional 22.32 cfs of ground-water discharge may 
bewacded sto. the previous. totalso£ -21./5eefs. j2thesesestimates 
Suggest that a grand total of 44.07 cfs of ground water™s 
discharging into Lake Ontario from the eastern hydrogeologic 
region. 
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SUMMARY 


A hydrogeological evaluation of the IFYGL Moira River, 
Wilton Creek and Thousand Islands study areas was made for the 
IFYGL porgram in order to estimate the amount of ground-water 
inflow to Lake Ontario. These study areas are considered to be 
representative of the hydrogeologic region extending from the 
Trent River to the eastern boundary of the Lake Ontario drainage 
basin at the St. Lawrence River near the Town of Prescott 
(Figure 1). The hydrogeology of the study areas was assessed 
utilizing field investigations of the geology, data from obser- 
vation-well networks and information from water-well records on 
file with the Ministry of the Environment. 

The IFYGL Moira River lies in parts of the physiographic 
regions of the Napanee Plain and Prince Edward Peninsula (Map 1). 
The physiographic regions of the Napanee Plain and the Leeds 
Knobs and Flats wholly enclose the IFYGL Wilton Creek and 
Thousand Islands study areas, respectively (Map 1). Regional 
topographic gradients are in the order of 33 feet per mile 
(0.0065 feet per foot) in the IFYGL Moira River and Wilton Creek 
study areas and in the order of 40 feet per mile (0.0075 feet 
pert foot) in the EFYGE Thousand Islands#study area. 

The climates in the study areas are influenced by the 
presence of Lake Ontario and as a result are very similar. Long- 
term, annual meteorologic data (tables 1, 2) suggest that the 
climate of the hydrogeologic region is slightly wetter toward the 
east and that the temperature is relatively uniform in all the 
areas under consideration. 

The bedrock in the hydrogeologic region consists of 
crystalline rocks of Precambrian age overlain by clastic rocks 
of Cambrian and Middle Ordovician age (Potsdam Formation and 
Shadow Lake Formation, respectively) in turn overlain by marine 
carbonates of Middle Ordovician age (Table 3). A veneer of 
glacial debris is found capping the bedrock over parts of the 
area but ror the most part bedrock is found outcropping - sous 
surface information from water-well records on file with the 
Ministry of the Environment indicate that the configuration of 
the bedrock surface is similar to that of the land surface 
(Map 2). The drainage initiated on the bedrock surface appears 
to be coincident with the present-day drainage and seems to be 
Controlled, for the most part, by structural features. (joimcem 
Eractures, etc.) in the bedrock. The bedrock 1s the primary, 
source of water for domestic and agricultural uses in the study 
areas under consideration. 

An overall decrease in the storage of the ground-water 
reservoir occurred in the study areas during the Field Year 
pertod from April 1, 1972 to March 31, 1973. “This is*imdteatea 
by an average decrease in water levels of 0.97, 0.18 and 0.28 
feet (Table 4) in the observation-well networks for the IFYGL 
Moira River, Wilton Creek and Thousand Islands study areas, 
respectively, Climatological data (Table 5) indicate that excess 
precipitation of 3.67, 6.47 and 5.23 inches over the long-term 
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annual average occurred in the IFYGL Moira River, Wilton Creek 
and Thousand Islands study areas, respectively during the Field 
Year period, These data suggest that a decrease in precipitation 
over the long-term annual average will probably result ina 
shortage of ground-water supplies in the study areas, 

The variability in the magnitude of the observation-well 
hydrographs from the observation-well networks in the study 
areas (Figure 3) suggest that a significant variation in the 
hydraulic properties of the bedrock is present, 

Water-level data from water wells completed in different 
water-bearing horizons in different seasons and years suggest 
that the general direction of ground-water flow can be determined 
from these data (Map 5). Ground-water divides are generally 
coincident with surface-water divides indicating that ground- 
water movement is from topographically high to topographically 
low areas, 

Well-productivity data suggest that the yield of bedrock 
wells in the study areas is independent of the lithology and 
primarily dependent on the fracture density of the bedrock, 
Wells completed in either the marine carbonates and claystone of 
Middle Ordovician age or crystalline rocks of Precambrian age 
Mendeminesnecirtic capacity from 00014 to 30 gpm per toot.ot 
drawdown with an average specific capacity of approximately 
0.20 Igpm per foot of drawdown (Figure 4), 

In order to evaluate the quantity of ground water discharging 
to Lake Ontario, an estimate of the hydraulic conductivity or 
transmissivity of the materials through which the ground water 
flows is required, Hydraulic conductivities of the overburden 
materials and bedrock were computed using the modified, non- 
equilibrium well formula and also using specific-capacity data 
(Table 7), Hydraulic conductivities derived from specific- 
Capacity data are larger by a factor of four to six than those 
hydraulic conductivities derived from the modified, non- 
equilibrium well formula, 

The contributing thickness of bedrock to significant ground- 
water flow was estimated from water-well data to be 50 feet in 
the IFYGL Moira River and Wilton Creek study areas, and 100 feet 
in the IFYGL Thousand Islands study area, The latter estimate 
of 100 feet suggests that fractures extend to a greater depth in 
the harder crystalline rocks of the IFYGL Thousand Islands 
study area than those fractures developed in the carbonate 
terrain of the IFYGL Moira River and Wilton Creek study areas, 

Using hydraulic conductivity values derived from short- 
term pumping tests, an average transmissivity of about 290, 
1500 and 3400 Igpd per foot was computed for the materials near 
the Lake Ontario shore in the IFYGL Moira River, Wilton Creek 
and Thousand Islands study areas. Ground-water inflow to Lake 
Ontario was estimated to be 0.155, 0.093 and 0,249 cts per mite 
in these study areas, respectively, giving a total of 
approximately 22 cfs for the hydrogeologic region that extends 
eastward from the Trent River to the eastern boundary of the 
Lake Ontario drainage basin on the St. Lawrence River. An 
additional amount of direct, ground-water discharge to Lake 
Ontario was computed for the shoreline of Price Edward County 
giving a grand total of approximately 44 cfs of ground-water 
discharge from this hydrogeologic region, 
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Map 1-Physiographic Regions, Meteorological Stations 
and Observation Wells 


Map 2 - Bedrock Geology and Topography 
Map 3-Locations of Water Wells 

Map 4 - Surficial Geology 

Map 5 - Generalized Potentiometric Surface 
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